Mon. Not. R. Astron. Soc. 000, CHE] (2006) 



Printed 5 February 2008 



(MN M£X style file v2.2) 



An improved model for the formation times of dark matter 
haloes 

Carlo Giocoli 1 ' 2 , Jorge Moreno 1 , Ravi K. Sheth 1 & Giuseppe Tormen 2 * 

1 University of Pennsylvania, 209 South 33 rd Street, Philadelphia, PA 19104 - 6396, USA 

2 Dipartimento di Astronomia, Universita degli Studi di Padova, Vicolo dell'osservatorio 2 1-35122 Padova, Italy 



ABSTRACT 

A dark matter halo is said to have formed when at least half its mass hass been 
assembled into a single progenitor. With this definition, it is possible to derive a simple 
but useful analytic estimate of the distribution of halo formation times. The standard 
estimate of this distribution depends on the shape of the conditional mass function — 
the distribution of progenitor masses of a halo as a function of time. If the spherical 
collapse model is used to estimate the progenitor mass function, then the formation 
times one infers systematically underestimate those seen in numerical simulations of 
hierarchical gravitational clustering. We provide estimates of halo formation which 
may be related to an ellipsoidal collapse model. These estimates provide a substantially 
better description of the simulations. We also provide an alternative derivation of the 
formation time distribution which is based on the assumption that haloes increase 
their mass through binary mergers only. Our results are useful for models which relate 
halo structure to halo formation. 
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1 INTRODUCTION 



Large samples of galaxy clusters selec ted in the optical 
iMiller et al.l 120051 ; iMcKav et all 120051 ) and other bands 
(ACT, SPT) will soon be available. These cluster catalogs 
will be used to constrain cosmological parameters. The tight- 
ness of such constraints depends on the accuracy with which 
the masses of the clusters can be determined from observed 
properties. These properties are expected to depend on the 
formation histories of the clusters. Models of cluster for- 
mation identify clusters with massive dark matter halos, so 
understanding cluster formation requires an understanding 
of dark halo formation. 

There is also some inte rest in using the d istribution of 
galaxy velocity dispersion s llSheth et al.ll2003l) to constrain 
cosmological parameters (N ewman fc David II20011 ). The ve- 
locity dispersion of a galaxy is expected to be related to 
the concentration of the halo which surrounds it, and this 
concentration is expecte d to be influenced by the forma- 
tion history of the h alo (|Tormenlll99Sl ; iBullock et alll200ll ; 
IWechsler et al.ll2002l ). Hence, this program also benefits from 
understanding the formation histories of dark matter halos. 

In hierarchical models, the formation histories of dark 
matter halos are expected to depend strongly on halo mass — 
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massive halos are expec ted to have formed more recently 
l|Press fc Schechteilll974l ). But quantifying this tendency re- 
quires a more pre cise definition of wha t one means by the 
'formation time '. lLacev fc Coll l|l993l ) provided a simple 
definition — it is the earliest time when a single progenitor 
halo contains half the final mass. For this definition, they 
showed how to estimate the d istribution of halo formation 
times. ISheth fc Tormenl \ 2004) provide associated expres- 
sions for the joint distribution of formation time and the 
mass at formation. This estimate depends on the distri- 
bution of progenitor masses at earlier times, and Lacey & 
Cole used the assumption that halos form from a spher- 
ical collapse to estimate this conditional mass function. 
However, a model based on ellipsoidal collapse provides a 
mor e accurate description of th e abundances of dark ha- 
los llSheth. Mo fc Tormenl 1200 lh and of their progenitors 
jSheth fc Tormenll2002l )." Hence, one expects to find that the 
ellipsoidal collapse model also provides a better description 
of halo formation. The primary goal of the present work is 
to study if this is indeed the case. 

Section [5] reviews Lacey & Cole's argument, and de- 
scribes how we estimate the ellipsoidal collapse based pro- 
genitor mass functions. Explicit formulae are provided in 
Appendix fX] Section [3] compares the model predictions for 
halo formation with measurements in simulations of hier- 
archical gravitational clustering. The simulations, and how 
we estimate halo formation times in them, are described in 
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Section r3.ll Sections 13 . 21 and 13 . 31 show that the spherical col- 
lapse model predicts lower formation redshifts than are seen 
in the simulations, particularly at low masses. Insertion of 
ellipsoidal collapse based expressions in Lacey & Cole's for- 
malism provides substantially increased accuracy. However, 
the agreement is not perfect — although the ellipsoidal col- 
lapse based expressions provide a good description of how 
the median formation redshift decreases with halo mass, it 
predicts formation time distributions which are broader than 
seen in the simulations. A final section summarizes our re- 
sults and discusses possible reasons for the discrepancy. 

There is some confusion in the literature about the re- 
lation between halo formation as defined by lLacev fc Cold 
l| 19931 ). and a quantity which arises in binary merger models 
of hierarchical clustering. The relation between this other 
quantity and the one defined by Lacey & Cole is clarified in 
Appendix [Bl 



2 THE FORMATION TIME DISTRIBUTION 

Following lLacev fc Cole! l| 19931 ). we will define the formation 
time of a halo as the earliest time that at least half of its 
mass has been assembled into a single progenitor. 

2.1 Relation to the progenitor mass function 

Consider an ensemble of halos of mass M at time T, and let 
N(m,t\M, T) denote the average number of progenitors of 
these halos that have mass m < M at time t < T. Since a 
halo can have at most one progenitor of mass m > M/2, the 
fraction of halos which have a progenitor of mass m > M/2 
at time t is 



/ < 

J M/2 



dmN(m,t\M, T). 



But, because they have a progenitor of mass m > M/2, 
these halos are also the ones which formed at some U < t. 
Hence, 



Jo 



dtp(t\M,T) 



dmN(m,tf\M, T), 
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where p(t\M, T) denotes the probability that a halo of mass 
M at T formed at time t. Differentiating with respect to tf 
yields 

d (' M 

p(t t \M,T) = — / dmN(m,t t \M,T) (2) 
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(3) 



Evidently, the formation time distribution is closely related 
to the distribution of progenitor masses and its evolution. 
Different estimates of the progenitor mass function will re- 
sult in different formation time distributions. 

In what follows, we will estimate the distribution of halo 
formation times using equation ([2]) . Thinking along the lines 
of equation instead shows how p(tf\M, T) can be related 
to quantities which arise naturally in binary merger models 
of hie rarchical clus t ering . Appendix iBl provides details. 

lLacev fc Coiel ()l993l ) used estimates for N(m, t\M, T) 
which were derived from a model in which halos form from 



a sph erical collapse l|Gunn fc Gottlll972l ; iPress fc Schechterl 
1974). In this case, they showed that the distribution of 
formation times could be scaled to a universal form: 



p(u>) = 2ui erfc(w/v / 2) 



where 



S B c(zf) — S sc (z ) 
s/S{M/2) - S{M) 



(4) 



(5) 



and 5 = 1. Here S sc (z) is the overdensity required for spher- 
ical collapse at z, and S(M) is the variance in the linear 
fluctuation field when smoothed with a top-hat filter of scale 
R — (3Af/47rp) 1//3 , where p is the comoving density of the 
background. In essence, u) is simply a scaled time variable: in 
an Einstein de-Sitter background cosmology us oc (zf — zq), 
where the constant of proportionality depends on the final 
mass M. 

Strictly speaking, this expression is valid for a white- 
noise power spectrum (S(m) oc m _1 ), but it has been found 
to provide a reasonable aproximation for more general power 
spectra as well. For this distribution, the median value of u> 
is 0.974, and the associated cumulative distribution is 



P{><j) 



2 /2 



+ (1-^ 2 ) 



2^ t I w 
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(6) 



However, recent work has shown that the spherical col- 
lapse model predicts fewer massive halos and more interme- 
diate mass halos than ar e seen in simulations of hierarchical 
gravitational clustering l|Sheth fc Tormenlll999h . Models in 
which hal os form from an ellipsoidal c ollapse may be more 
accurate (|Sheth. Mo fc Tormenl l200ll ) . These models also 
provide a better descrip tion of the progenitor mass function 
ijSheth fc Tormenl 120021 ). In what follows, we will show the 
effects of substituting the ellipsoidal collapse based expres- 
sions for N(m,t\M,T) in equation ([1]). We will also show 
that simply taking values of q smaller than one will make 
equation @ a good fit to the formation redshift associated 
to the ellipsoidal collapse. 

2.2 Progenitor mass functions in the excursion 
set approach 

We will use two approximations for the progenitor mass 
function. Both forms are derived from casting the ellip- 
soidal collapse model in the same language u sed for the 
spheri cal model — the excursion set formalism of lBond et al.l 
(| 1991T ) . In this formalism, halo abundances at a given time 
are associated with the first crossing distribution by Brow- 
nian motion random walks, of a barrier whose height de- 
creases with time, and may in addition depend on how many 
steps the walk has taken. Figure [1] illustrates: it shows a 
walk which starts at some initial S(M), and walks to the 
right, eventually crossing the barrier at some larger value 
of 5(771) > S(AI). Here S(m) is the variance in the ini- 
tial fluctuation field when smoothed with a tophat filter of 
comoving scale R = (3m/47rp) 1//3 , (p is the comoving back- 
ground density), and the height of the walk represents the 
value of the initial density field when smoothed on scale R. 
In hierarchical models, S(m) is a monotonically decreasing 
function of mass m. Hereafter we will use So and Si to in- 
dicate the overdensity associated with the spherical collapse 
model respectively at redshift zq and z\. 
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Figure 1. An example of a random walk and the three barrier 
shapes we consider in this paper. Here So refers to the critical 
value associated with spherical collapse overdenisity at redshift 
20. The jagged line is a sample brownian walk absorbed by the 
barrier associated with ellipsoidal collapse (solid curve). Short- 
dashed curves show the square-root barriers which we use to ap- 
proximate the ellipsoidal collapse barrier, and horizontal dotted 
lines show the constant barrier associated with spherical collapse. 
The upper set of barriers are associated with collapse at an earlier 
time. 



In the excursion set approach, the shape of the barrier, 
B(S), depends on the collapse model. If the barrier is crossed 
on scale S(m), then this indicates that the mass element 
on which the sphere is centered will be part of a collapsed 
object of mass m. Bond et al. used the fraction of walks 
which cross B(S) at S(m) as an estimate of the fraction 
of mass in halos of mass m: the parent halo mass function. 
Similarly, the fraction f(s\S) of walks which start from some 
scale S(M) and height B(S) and first cross the barrier B(s) 
at some s > S, can be used to provide an estimate of the 
progenitor mass function of M halos: 



.\tm\M Kim = ( — ) /(s|>'i.U- 



(7) 



(note that m < M because s > S). Thus, in the excursion set 
approach, the shape of N(m\M) depends on how the shape 
of the barrier changes with mass — it is in this way that the 
collapse model affects the parent halo and progenitor mass 
functions. 

The barrier shape is particularly simple for the spheri- 
cal model: B = S sc is the sa me constant for all values of S. 
ISheth. Mo fc Tormenl l|200lf ) argued that the barrier associ- 
ated with ellipsoidal collapse scales approximately as 

~S(m) 



B(S) = y/qSn 11 + 13 



q5l 



(8) 



where S sc is the value associated with complete collapse in 
the spherical model. Following Sheth et al., we set /3 = 0.5, 
7 = 0.6 and q — 0.707; the values of the first two parameters 



are motivated by an analysis of the collapse of homogeneous 
ellipsoids, whereas the value of q comes from requiring that 
the predicted halo abundances match those seen in simula- 
tions. The standard spherical model has /3 = and q = 1; 
we will also show the predictions of this model if (5 — but 
q = 0.707. 

For general 7, exact analytic solutions to the first cross- 
ing distribution of barriers given by equation © are not 
available. For this reason, we have studied two analytic ap- 
proximations which result from addressing the first crossing 
problem in two different ways. The first approach uses an an- 
alytic approximation to the first crossing distribution which 
ISheth fc Tormenl (|2002l) showed was reasonably accurate. 

Our second approach is to substitute the ellipsoidal col- 
lapse barrier for one which is similar, but for which an ex- 
act analytic expression for the first crossing distribution is 
available. Specifically, when 7 = 0.5, then the barrier height 
increases with the square-root of S, and the first crossing dis- 
tribution can be written as a sum of parabolic cylinder func- 
tions (Breiman 1966). In this approach, we approximate the 
progenitor mass function using an expression which would 
be exact for a square-root barrier, bearing in mind that the 
square-root barrier is an approximation to the ellipsoidal col- 
lapse dynamics. For this barrier shape, we set (3 = 0.5 and 
q — 0.55 since these choices result in parent halo and progen- 
itor mass functions which best fit the simulations; Figure [T] 
illustrates the similarity of the barrier shapes. We will show 
that both these approximations provide substantially bet- 
ter descriptions of the simulations than does the expression 
which is based on spherical collapse. Explicit expressions 
for these barrier crossing distributions are provided in Ap- 
pendix^] They are related to the progenitor mass function 
by equation (0 . 



3 COMPARISON WITH SIMULATIONS 

The simulation data analysed below are publicly available 
at http://www.mpa-garching.mpg.de/Virgo, The simula- 
tio n itself is known as GIF2, and is described in some detail 
bv lGao et ail (|2004l ). It is of a flat ACDM universe with pa- 
rameters (n m ,vg,h,n b h 2 ) = (0.3,0.9,0.7,0.0196). The ini- 
tial fluctuation spectrum had an inde x n = 1, with trans- 
fer fu nction produced by CMBFAST ()Seliak fc Zaldarriaeal 
Il996j) . 



3.1 The GIF2 simulation 

The GIF2 simulation followed the evolution of 400 3 particles 
in a periodic cube 110 /i _1 Mpc on a side. The individual par- 
ticle mass is 1.73 x 10 9 /i -1 M@, allowing us to study the for- 
mation histories of haloes which are one order of magnitude 
smaller in mass than was previously possible. Particle posi- 
tions and velocities were stored at 50 output times, logarith- 
mically spaced between 1 + z = 20 and 1 + z — 1. Halo and 
subhalo mer ger trees were c o nstru cted from these outputs as 
described bv lTormen et all (|2004l ). Briefly, halos at a given 
redshift are identified as spherical objects which enclose the 
virial density appropriate the cosmological model at that 
redshift (computed using the spherical collapse model). 

A merger history tree was constructed for each halo 
more massive than 10 11 ' 5 Ii~ 1 Mq (i.e., containing more than 
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Figure 2. Cumulative distribution of dark halo formation times for halos identified at z = 0. From top to bottom, panels show results 
for halos with masses in the range logioAf//i _1 Mq : 11.5-12, 12-12.5, 12.5-13, 13-13.5, and 13.5-14. Symbols show the measurements in 
GIF2; dotted curve shows the prediction associated the constant barrier spherical collapse model; dot-dashed curve shows the analytical 
fit equation (0 with q = 0.707. Short-dashed and solid curves show the predictions associated with the square-root and ellipsoidal 
collapse based models. 



180 particles) as follows. The progenitors of a halo of mass 
Mo at zo are found by identifying in the previous output 
zi > zo all haloes which contribute any number of particles 
to Mo- Of these, the progenitor which contributes the most 
mass to Mo is called Mi. This procedure is then repeated, 
starting with Mi at z\ , considering its progenitors in the pre- 
vious output time 22 > z\ , and choosing from them the one, 
Mi , which contributes the most mass to Mi . In this way, the 



mass of the most massive progenitor is traced backwards in 
time, to high redshift. 

We will present our results as a function of halo mass. 
We have estimated formation times for 5611 halos with 
logi O (M//i _1 M0) in the range 11.5 - 12, 2431 in the range 

12 - 12.5, 892 in the range 12.5 - 13, 341 in the range 

13 — 13.5, 92 in the range 13.5 — 14, and 29 in the range 

14 — 14.5. These abundances are well fit by the formula of 
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Figure 3. Same as previous figure, but now for halos identified at z = 0.5. As before, the spherical model predicts lower formation 
redshifts than observed. Whereas the ellipsoidal collapse based expressions predict the median formation rcdshift reasonably well, they 
predict a broader range of redshifts than are observed. 



ISheth fc Torn^nl (|l999l ). so they are in better agreement 
with the ellipsoidal rather than spherical collapse model. 



3.2 Cumulative distribution of formation times 

Figure [2] shows the cumulative distribution of dark halo 
formation redshifts (i.e., equation [TJ for halos identified at 
z = 0. In general the analytical expressions obtained from 
the two different approxima te solutions to the ellips oidal col- 
lapse barrier problem — the (|Sheth fc Tormer]|2002l ) approx- 



imation for the first crossing distribution of the ellipsoidal 
collapse barrier, or the exact expression for the first cross- 
ing distribution of the square root barrier — are complicated. 
However, we have found that the predictions of the two el- 
lipsoidal based models are quite well approximated by the 
expression for the spherical model with n — 0, equation ©, 
by simply changing the value of q. 

Different panels show results for the mass bins described 
above. The points show measurements in the GIF2 simula- 
tion, and the four curves show the formation time distri- 
butions associated with the ACDM spherical collapse (dot- 
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Figure 4. Median formation redshift for halos identified at z = 
as a function of halo mass. Points with (Poisson) error bars 
show our measurements in the GIF2 simulation. Smooth curves 
show the median formation times associated with three different 
models of halo formation: spherical collapse (dotted), ellipsoidal 
collapse (solid) and the square-root barrier approximation (short- 
dashed). Dot-dashed curve shows the prediction of equation IllOjl 
with q = 0.707. 

ted), ellipsoidal collapse (solid), square-root barrier approx- 
imation (dashed) and with equation © using q — 0.707 
(dot-dashed). 

Again, the spherical collapse model severely underesti- 
mates the redshifts of halo formation. The two ellipsoidal 
collapse based estimates fare better, in the sense that they 
predict median formation redshifts which are closer to those 
seen in the simulation. 

However, both ellipsoidal based estimates clearly pre- 
dict a wider range of formation redshifts than is seen in the 
simulation — at fixed mass, the distribution of halo formation 
redshifts is narrower than predicted. This remains true for 
the formation time distribution of halos identified at z — 0.5 
shown in Figure [3] 

3.3 Mass-dependence of median formation time 

Figure [4] shows the median formation redshift of halos iden- 
tified at z = as a function of halo mass. Points show 
our measurements in the GIF2 simulation The simulation 
showes that massive halos formed more recently: the me- 
dian formation redshift decreases with halo mass. 

Smooth curves show the formation time distribu- 
tions associated with three different models of halo forma- 
tion: ACDM spherical collapse (dotted), ellipsoidal collapse 
(solid) and the square-root approximation (short-dashed). 
In all cases, the model predictions for the median formation 
redshift were obtained by finding that Zf at which 




Figure 5. Same as previous Figure, but for halos identified at 
z = 0.5 in the GIF2 simulation. 

The Figure shows that halos less massive than 10 14 /i _1 Mq 
clearly form at higher redshifts than predicted by the spher- 
ical collapse model (dotted line). Both our estimates of 
the ellipsoidal collapse prediction are in substantially better 
agreement with the simulations. The prediction associated 
to equation ([6]) is given implicitely by inserting the median 
rescaled formation redshift Q = 0.974 into equation ((SJ), pro- 
viding: 

974 / 

5 ac (zi\M,z ) = S BC {z ) + — — y/S(M/2) ~ S(M). (10) 

The choice q — 0.707 is shown by the dot-dashed line in the 
Figures; setting q = 0.6 yields even better agreement with 
the simulations. 

Fitting functions, accurate to a few percent, for 5 sc {z) 
and S(m) are available in the literature. For instance, 

s BC (z) « d+(z) A ( i2^) 2 / 3 n° m 0055 (11) 

l)Navarro et alJ fl997t ) . where Q, m is the ratio of the back- 
ground to the critical density at z, and 

D+[z) = ^fi m [^ /7 -n A +(l+in m ) (i+ JLfi A )] ~* (12) 

jCarroll et al.lll992T ). and, for a CDM power spectrum, 

S(m) = A (1 + 2.208m p - 0.7668m 2p + 0.7949m 3p ) " 4/(9p) 

(13) 

ijTaruva fc Sutoll2000l ). where A is the normalization factor 
of the linear theory power spectrum at z = (so it depends 
on as), fh = m(hT) 2 /10 12 M@, where F is the parameter 
which describes the shape of the power spectrum (typically 
rwfift), and p = 0.0873. 



4 DISCUSSION AND CONCLUSIONS 

lLacev fc Cole] l| 19931 ) defined the formation time of an object 
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as the earliest time when at least half its mass was assembled 
into a single progenitor. We clarified the relation between 
this definition of formation time, and a quantity which arises 
in binary merger models of clustering (Appendix |B} . We 
have shown that insertion of spherical collapse based expres- 
sions in Lacey & Cole's (1993) formalism for halo formation 
underestimates the redshifts of halo formation; this is consis- 
tent with previous work (Lin, Jing & Lin 2003). Ellipsoidal 
collapse based expressions are a marked improvement: al- 
though they result in formation redshift distributions which 
are broader than seen in simulations (Figures [3] and©, they 
predict the median formation redshift quite well (Figure [4] 
and©. 

The fact that our predicted formation time distribu- 
tions are broader than those seen in the simulation can 
be traced back to the fact that the low redshift progeni- 
tor mass functions from the excursion set approach are not 
in particularly good agreem ent with the simu l ations . This 
is similar to the findings of ISheth fc Tormenl (2002), who 
noted that some of the discrepancy was almost certainly 
due to the idealization that the steps in the excursion set 
walks are uncorrelated. Recent work has shown that there is 
some correlation between halo formation and environment 
dSheth fc Tormenll2004l : iGao et al.ll2005l ; lHarker et al"1l2006l : 
IWechsler et al.l[200rt ) — this almost certainly indicates that 
a model with uncorrelated steps will be unable to provide a 
better description of the simulations. Nevertheless, the fact 
that the median formation redshifts are quite well repro- 
duced by our model does represent progress. 

Although the exact expressions associated with our 
models are complicated, we were able to find a useful fitting 
formula for the median formation redshift, equation (I10p , 
which we hope will be useful in studies which relate the for- 
mation times of halos to observable quantities. 
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APPENDIX A: EXPLICIT EXPRESSIONS FOR 
BARRIER CROSSING DISTRIBUTIONS 

We are interested in walks which start from one barrier 
Bq(S) and walk to another, -Bi(s), with s > S. Thus, the 
barrier to be crossed is of the form 

B 1 (s)-B (S) = [l+/?(-^) 7 ] -^S [l+/?(Jg)^ ]. 

°(A1) 

In the formulae do and S\ refer to the sferical collapse over- 
denisity at redsh i ft zp and z\. For barriers of this form, 
ISheth fc Tormenl (|20021 ) showed that 

ds \T(s\S)\ 



f(s\S)ds « 
where 



8 ^2tt(s - S) 

[Bi(s) - flo(S)] 2 



" XP > 2(s-S) 

(S - 8) n d n [B!(8) - Bp(S)] 

n! ds n 



(A2) 



(A3) 
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When ,3 = and q = 1, t his reduces to the sph erical collapse 
based expression used bv lLacev fc Cole] (|l993l ): 



f(s\S)ds 



ds 



Si — 5 



exp 



(Si - So) 2 
' 2(s-S) 



(A4) 



- S ^/2tt(s - S) 
When 7 = 1/2, then 

Bi(s) - B (S) = s/q{5i - So) - /3Vs + f3^, (A5) 
and the exact solution to this square-root barrier problem is 



f(s\S)ds 



exp 



E 



4 - P 2 



D x (w p )D' x (-0) ds 



Ix{-/3) 2s 

(|Breimanlll96rj ; lMahmood fc Raieshll2005l ), where 
^/q(5i - So 



(A6) 



Wf) 



dx D x (a 



Vs j-0 

and D\(x) axe parabolic cylinder functions which satisfy 

D' x (x) = (x/2) D x (x)-D x+ i(x). 
The sum is over the roots A defined by 

£>*(-/?) = 0, so D' x (-l3) = -D x+ i(-f3). 



APPENDIX B: ON THE DIFFERENCE 
BETWEEN HALO FORMATION AND 
CREATION 

The main text dea ls with halo formation as defined by 
lLacev fc Coll (|l993h . Before continuing, it is worth men- 
tioning that there is another context in which the term 'halo 
formation' arises — this is when the halo population is mod- 
eled as arising fro m a binary coagulatio n process of the type 
first described by I SmoluchowsklPl 1 9 1 (i ) . In this description, 
the time derivative of the halo mass function is thought of 
as the difference of two terms: one represents an increase in 
the number of halos of mass m from the merger of two less 
massive objects, and the other is the decrease in the num- 
ber of m-haloes which results as m-haloes themselves merge 
with other halos, creating more massive halos as a result: 



dn(m, t) 
dl 



= C(m,t) - D(m,i) 



(Bl) 



The gain term, the first on the rig ht hand side above, is 
sometimes c alled the halo formation ijSheth fc~ PitmarJll997h 
or creation (jPercival fc Miller! 1 19991 : IShethl l2003h term. In 
what follows, we will use the word 'creation' to mean this 
term, and 'formation' to mean the quantity studied by Lacey 
& Cole. 

Halo formation and creation are very different quan- 
tities, as we show below. Nevertheless, the y are some- 
times used interchangeably in the literature fe.g. lVerde et al.l 
(2001)). Here we show explicitly how to compute the Lacey- 
Cole formation time distribution from the Smoluchowski- 
type creation and destruction terms, with the primary aim 
of insuring that the error of confusing one for the other is 
not repeated. 

In the Lacey-Cole picture, the formation time distri- 
bution of A/-halos identified at T is given by equation ©. 



In this case, one first integrates over the conditional mass 
function, and then takes a time derivative. Here, we will in- 
stead take the time derivative inside the integral over mass 
and compute it before integrating over mass, as suggested 
by equation ©. In this case, the integrand has the form of 
a time derivative which plays a central role in the Smolu- 
chowski picture. In particular, we can write 



dN{m,t\M,T) 
dt 



= C(m, t\M, T) - D(m, t\M, T), (B2) 



where C and D are now the creation and destruction terms 
associated with the p rogenitor mass function, evaluated ex- 
plicitly in lShethl (|2003l ). 

This simple step shows explicitly that p(t\M, T) can be 
written in terms of Smoluchowski-like quantities: 



p(t\M,T) 



dmC(m,t\M, T)- 



M/2 



M/2 



dmD(m,t\M, T). 

(B3) 

A little thought shows why this works out so easily. Sim- 
ply integrating the halo creation rate C(m,t\M, T) over the 
range M/2 ^ m ^ M overestimates halo formation, since 
some of the halos created at t, with mass 3M/4 m M 
say, may actually have been may created by binary merg- 
ers in which one of the halos had a mass in the range 
M/2 ^ m ^ 3M/4. Such creations should not be counted 
towards halo formation, since formation refers to the ear- 
liest time that m exceeds M/2. However, it is precisely 
this double-counting which the second term, the integral of 
D(m, t\M, T) over the range M/2 ^ m ^ M, removes. Note 
in particular that, whereas the formation time distribution 
p(t\M, T) is related to Smoluchowski-type quantities, it is 
not the same as C(m,t). For completeness, note that 



C{m,t) 



f 

■J 111 



dMn(M, T)C(m,t\M,T), 



(B4) 



where n(M, T) denotes the number density of halos of mass 
M at time T. 
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